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The motion of an initially spherical capsule in a wall-bounded oscillating shear flow is 
investigated via an accelerated boundary integral implementation. The neo-Hookean 
model is used as the constitutive law of the capsule membrane. The maximum 
wall-normal migration is observed when the oscillation period of the imposed shear 
is of the order of the relaxation time of the elastic membrane; hence, the optimal 
capillary number scales with the inverse of the oscillation frequency and the ratio 
agrees well with the theoretical prediction in the limit of high-frequency oscillation. 
The migration velocity decreases monotonically with the frequency of the applied 
shear and the capsule-wall distance. We report a significant correlation between 
the capsule lateral migration and the normal stress difference induced in the flow. 
The periodic variation of the capsule deformation is roughly in phase with that of 
the migration velocity and normal stress difference, with twice the frequency of the 
imposed shear. The maximum deformation increases linearly with the membrane 
elasticity before reaching a plateau at higher capillary numbers when the deformation 
is limited by the time over which shear is applied in the same direction and not by 
the membrane deformability. The maximum membrane deformation scales as the 
distance to the wall to the power 1/3 as observed for capsules and droplets in near¬ 
wall steady shear flows. 
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I. INTRODUCTION 


The dynamics of a capsule in external flows have attracted enormous interest due to their 
physiological and biological significance, leading to the discovery of a variety of interesting 
and complex phenomena. Since the observation of the tank-treading (TT) motion of a 
human red blood cell (RBC) in the viscometer ——, real cells and their models have been 
examined in shear flow, one of the simplest flows, both experimentally, theoretically and 
numerically in order to understand the motions of biological cells^r—. The focus was on 
an initially spherical capsule that undergoes shear-induced TT motion, reaching a steady 
ellipsoidal shap e 17,18 , when the membrane elements rotate along the stationary configuration 
like the caterpillar driving of a tank. The undeformed shape of biological cells is, however, 
seldom a sphere, rather geometrically anisotropic, the bi-concave shape of undeformed RBCs 
being a classic example. This non-sphericity introduces a shape memory for the capsules and 
tends to orient them in a preferential direction—. As a consequence, the capsules tank-tread 
and oscillate around a preferred direction as the shear rate is above a certain value^ff^L; this 
oscillatory behavior is termed as the swinging mode^-. 


Oscillating shear flows have been investigated to account for the physiological pulsation, 
and more complex and diverse capsule dynamics explored. Experiments^ showed that the 
RBCs deform more in the retarding phase than in the accelerating phase. Cells in har¬ 
monically modulated shear were examined by theoretical approaches^ to reveal a resonance 
behavior; a particular combination of the oscillation frequency and phase can induce the 
tumbling motion of a capsule which would otherwise swing under steady shear flow. Dupire 
et ai— theoretically produced disordered motions of the RBCs under sinusoidally varying 
flows with physiologically relevant parameters. A similar chaotic behavior was also repro¬ 
duced in the analytical investigation by Noguchi ™ who found multiple limit-cycle solutions 
of a model RBC subject to an oscillating shear flow at high frequency; their tumbling or 
tank-treading motions were found at high or low shear amplitudes. The recent numerical 


simulations in Ref. 


271 reproduce the two swinging modes of an initially oblate capsule by 


varying the frequency of the external shear; these authors also document the high sensitiv¬ 
ity of the motion of the capsule to its initial orientation, indirectly confirming the chaotic 
motion reported in Ref. 


25 


These previous studies mainly focused on the dynamics of capsules in an unbounded os¬ 
cillatory shear flow. In the micro-circulation, however, cells like the RBCs move in strongly 
confined oscillating flows, an example being the physiological process called vasomotion. 
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This corresponds to the spontaneous and rhythmic oscillations of large arteries as well as 
microvessels22~— . Vasomotion is independent of the common physiological pulsations in¬ 
duced by the heart beat, innervation or respiration^, but due to the constriction and dila¬ 
tion of the smooth muscle. Its possible benefits include reducing the time-averaged hydraulic 
resistance to the blood flow22 and the consequent enhancement of oxygen transporfr^ 1 ^. It 
is therefore important to understand the dynamics of deformable particles like biological 
cells in a confined oscillating flow. 

In addition to the complexity introduced by the oscillation, the influence of wall confine¬ 
ment on the cell motion is crucial and has been, thus far, mostly investigated in steady flows 
where cells exhibit cross-streamline migrations in most cases. Fahraeus and Lindqvist— first 
unraveled the biological importance of confinement. They discovered that the RBCs tend to 
deform significantly and migrate towards the center of a microvessel owing to the hydrody¬ 
namic lift resulting from the interaction between the cellular deformation and the wall; this 
migration facilitates a lower flow resistance and a more effective mass transport. Simula¬ 
tions using the front-tracking method have been conducted by Doddi et ai— to investigate 
the cross-stream lift of a capsule towards the centerline in a Poiseuille flow; these authors 
observed that the migration velocity scales as the distance to the wall to the power 1/3 but 
varies non-linearly with the capillary number. The recent numerical study by Pranay et 
al.~ 7 , based on an accelerated boundary integral method as that used here, investigates the 
lateral migration of an individual and a suspension of initially spherical capsules in wall- 
bounded Newtonian and viscoelastic shear flows. These authors show that the migration 
of the capsule depends on the capillary number, and that it is attenuated by the addition 
of polymers. Migration has also been reported for vesicles: in an unbounded Poiseuille 
flow^ 1 ^, the migration velocity depends on the flow curvature, vesicle deformability and 
viscosity contrast of the fluid inside and outside the vesicle; in a wall-bounded shear, the 
migration velocity is proportional to the wall normal component of the particle stresslet and 
the presence of the wall delays its transition from the tank-treading motion to trembling 
and tumbling^. 

Recently, two groups^ ,4? have numerically examined the case of a near-wall single capsule 
transported in a Newtonian shear flow, using a front-tracking and a boundary integral 
method respectively. Singh et al— document the power-law relations between the migration 
velocity, capsule deformation and the capsule-wall distance, as previously done analytically 
for a droplet. These authors also propose a semi-analytical theory identifying two competing 
mechanisms influencing the lateral migration: the interfacial stresses and the viscosity ratio; 
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this theory agrees well with the simulation results for large capsule-wall distances. Nix 
et al.— pay more attention to the correlation between the deformation of a capsule and 
its migration. These authors find that an asymmetric capsule deformation reduces the 
migration velocity and this effect is compensated by an increase in the capsule stresslet, in 
turn enhancing its migration. 

Very recently, Matsunaga et al.— investigated the motion of capsule in an unbounded 
oscillating shear flow, observing an overshoot phenomenon when the maximum deformation 
of capsule can be larger than its deformation in a steady shear. In this manuscript, we extend 


the work in Ref. 


41 


, Ref. 


42 


and Ref. 


43 


by considering a near-wall individual capsule in 


oscillating plane shear flow. We hence introduce one additional factor, i.e., the unsteadiness 
of the background flow that is a typical feature of biological environments. We aim to 
understand the interplay between the wall and the flow unsteadiness, explaining its influence 
on the dynamics of a capsule. 

The paper is organized as follows. After the description of the problem setup, we give 
a brief introduction to our computational framework based on an accelerated boundary in¬ 
tegral method suited for general geometries. We firstly report the trajectories of capsules 
migrating from the wall and then investigate the dependence of the lateral migration velocity 
on the capsule capillary number, the capsule-wall distance and the frequency of the oscil¬ 
lation of the background shear. We further analyze the correlation between the migration 
velocity and the normal stress difference that arises from the viscoelastic effects due to the 
deformable particles. Finally, the deformation of the capsule is examined; the phase differ¬ 
ence between the deformation and the applied shear is measured to quantify the induced 
viscoelasticity. The paper ends with a summary of the main findings. 


II. PROBLEM SETUP AND NUMERICAL METHODS 
A. Problem setup 

We compute the motion and deformation of a capsule subject to a wall-bounded shear 
flow, see the sketch in figure [U The undeformed shape of the capsule is a sphere with radius 
a and center of mass initially located at a distance h^/a above the wall. A time-periodic 
harmonic shear flow 7 (t) = 7 max cos (cot) is applied. 7 max denotes the maximum value of 
the shear, and 00 is the oscillation frequency. The characteristic time of the flow is thus 
T = 2n/co. In the current work, we study a capsule as a model cell and reproduce its motion 
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i(t) = 7maxCOs(o;t ) 



FIG. 1. (Color online) We consider a capsule whose undeformed shape is a sphere of radius a, 
subject to an oscillating wall-bounded shear flow 7 ( t ) = 7 max cos (ut). The center of mass of the 
capsule is initially a distance h- m \ above the wall and its lateral migration, measured from the initial 
position, is denoted /ii at . The distance between the capsule and the wall is thus h cw = h - ia i + hi at . 


and deformation in the micro-circulation or in micro-fluidic devices. We therefore assume 
that the capsule is advected by a creeping flow and the Reynolds number Re, indicating the 
ratio of inertial over viscous forces, is zero. 

The capsule is considered as a fluid-filled droplet enclosed by an infinitely thin elastic 
membrane. The fluid inside and outside the capsule has the same density p and the same 
dynamic viscosity / 1 . For real biological cells, the viscosity of the fluid inside and outside is 
not necessarily the same; this is not taken into account here to limit the investigated pa¬ 
rameter space. The capsule is deformed by the fluid flow and the stress on its membrane is 
determined by the neo-Hookean constitutive law^. The membrane has a shear modulus G s 
and zero bending stiffness. Elastic stresses develop on the surface of the capsule due to its 
deformation; this stress modifies the surrounding flow in return. I 11 this fluid-structure inter¬ 
action problem, the viscous stress of the flow and clastic stress of the deformable membrane 
compete with each other, their ratio defined as the capillary number Ca = /i 7 ma X /G s . 
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B. Numerical method 


The membrane of the capsule is discretized into N Lagrangian points. Neglecting the 
inertia of the capsule, the elastic force balances the flow viscous force on its surface, 


P + (^out O’;,,) • n out 0, 


( 1 ) 


where p c is the elastic force per unit area exerted by the membrane, cr out , cx in are the stress 
tensor of the flow outside and inside the capsule, and n out is the unit outwards normal 
vector. Given the distribution of p e on the membrane, the flow obeys at each time the 
Stokes equations with N point forces exerted on the fluid: 


N 


Vp + pV 2 u = — f i§ i 


i =1 


V ■ u(x) = 0, 


( 2 ) 


where p and u are the pressure and velocity field, f \ denotes the force on the fluid at the 
position Xj and can be approximated by 


f, = 


PldSt, 


( 3 ) 


where S t represents the elemental patch around the fth Lagrangian point. We solve the 
governing flow equations with a boundary integral method, accelerated by the general ge¬ 
ometry Ewald method (GGEM) proposed by Hernandez-Ortiz et al.—. GGEM decomposes 
the Stokes solution into two parts: i) the short-ranged interactions computed by traditional 
boundary integral techniques; and ii) the long-ranged interactions handled by a mesh-based 
Stokes solver. In our implementation, we take the Stokes sub-solver of the open source 
software NEK5000^ as the mesh-based solver. We perform singular and nearly-singular 
integration for an accurate near-field solution, see Ref. [47|. A proper treatment of the two 
is crucial to achieve the necessary numerical accuracy.— 

The flow domain is infinite and bounded only by a plane wall. However, in the frame¬ 
work of GGEM, a computational domain of finite size is needed; here we choose 24 a(x) x 
24 a(y) x 15a(z). We impose periodic boundary direction in the streamwise (x) and spanwise 
(z) direction, while keeping the bottom wall stationary and imposing at the upper bound¬ 
ary a time-periodic velocity in the x direction. The domain is discretized by 4000 cubic 
spectral elements with 4x4x4 Gauss-Lobatto-Legendre points and the grid is refined in 
the region where the capsule moves. This strategy has been used to study a capsule in 
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an unbounded shear, producing results in agreement with published data obtained with a 
traditional boundary integral method^. 

Once the flow solution is available, the velocity is known at each Lagrangian point, u (xj). 
Due to the no-slip, non-penetrating boundary condition on the membrane of the capsule, 
the rate of change of the position of the ith point is given by 


dxi 

dt 


= u x,- 


(4) 


Given the new coordinates of the N Lagrangian points, we are able to calculate the elastic 
force per unit area on each point. The computation of the stress is based on the displacement 
with respect to the undeformed shape of the capsule, through the constitutive law of the 
membrane, the neo-Hookean model here. We use a global spectral method based on the 
spherical harmonics^ to represent the surface of the capsule and to solve for the elastic 
force p c . The advantage of this approach is twofold: i) high order spatial derivatives on the 
material points are computed with high accuracy, which is crucial for the calculation of the 
elastic force; ii) the same spectral discretization can be used for the boundary integration 
performed when solving the short-ranged hydrodynamic interactions. For the details of our 


implementation, the readers are referred to Refs. 


50 


51 


One of the important features is 


the capability of simulating deformable capsules in general geometries, as illustrated by the 
design of a deformability-based cell sorting device^ and a constricted micro-fluidic channel 
that can be used to infer the mechanical properties of cellular particles^. 


III. RESULTS 
A. Trajectories 

We firstly look at the trajectories of the capsules. The time evolution of the center of mass 
and surface centroid of the capsule are depicted in the shear plane (x—y plane) in figure [2] for 
three cases, {Ca,u/j max ) = (0.15,1), (Ca, cj/j max ) = (0.3,1) and (Ca, cj/j max ) = (0.3,2). 
Note that the center of mass is equivalent to the volume centroid, as opposed to the surface 
centroid, and the two centroids coincide for centersymmetric objects while they deviate as 
this symmetry is lost. In this study, the capsule symmetry is broken by the presence of the 
wall. The motion is a combination of oscillations in the streamwise (x) direction due to 
the oscillating background flow and wall-normal lateral migration due to the hydrodynamic 
lift. The trajectories of the two centroids, h\ at /a and h\ at display significant differences: the 
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x/a x/a 

(a) (b) 

FIG. 2. (Color online) The trajectories of (a) (x/a, hut/a) of the center of mass, and (b) 
[x/a, hi a t/a^ of the surface centroid on the shear plane, of the capsules (Ca,ui/ 7 max ) = (0.15,1) 
(solid), (Ca, o;/7 max ) = (0.3,1) (dashed) and (Ca,uj/ 7 max ) = (0.3,2) (dot-dashed). The capsules 
are released with an initial offset h- m \/a = 2. 

center of mass moves vertically in a reciprocal way while the surface centroid migrates away 
from the wall continuously. 

This is also illustrated in figure [3] where we display the time migration of the two centroids, 
hi at /a and hi at /a. They both vary periodically with the frequency of the imposed shear, 
but neither resembles a regular sinusoidal as the shear: the former increases with time non- 
monotonically, its migration velocity can be negative while the latter increases monotonically 
and hence the velocity is always positive. This difference reflects the asymmetry of the shape 
of the capsule. They nevertheless have approximately the same time-averaged value; after 
all, both centroids represent the motion of the capsule and thus coincide with each other in 
a statistical sense. In the following, we mostly study the time-averaged capsule migration 
and thus only report results for the center of mass. 

We further note that a solid spherical particle will not migrate in the wall-normal direction 
in Stokes flow, while its ellipsoidal counterpart will exhibit an oscillatory vertical drift, with 
a zero mean. Nonlinearity is the cause of the particle migration, either due to inertial 
effects, a viscoelastic fluid or to the deformable membrane as in our case. As an example 
of solid particles, we note that a sphere has been shown to migrate in a constant shear 
flow of a viscoelastic fluid^ due to the nonlinear relation between the fluid stress and the 
strain rate. Under oscillating shear—, the spherical particle follows a wagging trajectory, 
very similar to that we observe here. The physical picture is also analogous to that of 
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(a) 


(b) 


FIG. 3. (Color online) (a) Lateral migration /ii at /a, of the center of mass, and (b) /ii at /a, of the 
surface centroid of the same capsules in figure [2j The dotted curve indicates the imposed periodic 
shear 7 = 7 max cos (cot), arbitrarily scaled for visualization purposes. The shaded region indicates 
the time when the the instantaneous lateral migration /ii a t and the shear 7 reach the local minimum 
and maximum respectively. 


a microswimmer exploiting nonlinearity to achieve a net locomotion despite a reciprocal 
motion, an interesting case in the low-Reynolds-number swimming dynamics^. In the Stokes 
regime, but in a viscoelastic or non-Newtonian fluid, a reciprocal swimming pattern can 
lead to a net displacement^— , although the same strategy would not work in a Newtonian 

fiuid^— . 


We further notice that for the same oscillation frequency, (u/ 7 max = 1, the motion of the 
two capsules with Ca = 0.15 and Ca = 0.3 is similar despite the different deformability. 
Conversely, the trajectory of the same capsule, Ca = 0.3, varies significantly when the 
frequency of the imposed shear increases from oj/ 7 max = 1 to u>/ 7 max = 2. The streamwise 
domain spanned by the capsule reduces to about half when doubling u;/ 7 max , as expected 
since we keep the same maximum shear. The capsule lateral migration also diminishes 


considerably in this case, as shown in figure 3(b) 


B. Velocity of the capsule 

In this section, we analyze the migration and slip velocity of the capsule. The lateral 
migration velocity U\ at , is computed as Ui at = dh cw /dt and the slip velocity U s u p = 7 h cw — U x 
where h cw = hmi + ^iat is the distance between the capsule and the wall and U x the streamwise 
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U\ a t/ aj 



FIG. 4. (Color online) (a) Lateral migration velocity U\ a .t/a A f mgx and|(b)[ slip velocity U s y ip /a A f max . 


as a function of the dimensional time t/T for capsules with Ca = 0.0375 (solid curve), Ca = 0.15 
(dashed) and Ca = 0.6 (dot-dashed), released with an initial offset h- m \/a = 2; the oscillation 
frequency of the imposed shear uj/ 7 max = 5/3. The time-averaged migration and slip velocities 
are denoted by horizontal lines with circles (Ca = 0.0375), squares (Ca = 0.15) and diamonds 
(Ca = 0.6). The dotted curve indicates the applied oscillating shear 7 max cos (ut ), arbitrarily 
scaled for visualization purposes. 


velocity of the center of the capsule. The temporal evolution of f/i a t/«7ma X and U s n p /a A / n 


are depicted in figure 4(a) and figure 4(b) for three capsules with the same initial offset 
hini/a = 2 and frequency of the imposed shear, co/ A f max = 5/3. After an initial transient, 
the two velocities vary periodically in time. The transient is also observed in the case of 
steady shear—, before the capsule reaches a quasi-steady state; only then its deformation 
and velocity do not vary with the instantaneous capsule-wall distance h cw . In the case of 
periodic shear investigated here, the capsule motion reaches a quasi-periodic regime. Unless 
otherwise specified, we only examine its dynamics at this stage, when its net migration /i lat 
is negligible compared to the initial offset h ini and h ini ~ h cw . In the following, h ini will 
therefore be used to denote the capsule-wall distance. 


As clear from figure 4(b), the slip velocity and the background shear flow have different 
signs. This indicates that the capsule moves faster than the local background flow, an 
effect more pronounced when t/T approximately assumes integer values and the capsule 
returns about the initial location (x = 0). This is different to what is observed in a steady 
shear flow where the capsule always lags behind the surrounding flui d 41 ' 42 . In fact, as 
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FIG. 5. (Color online) Lateral migration velocity Li a t/ a 7 max as a function of the dimensional time 
t/T and the capillary number Ca, at fixed oscillation frequency ui/ 7 max = 5/3. Data are reported 
for three initial offsets, h-mi/a = 1.5, h-^/a = 2 and h\ n \/a = 3 indicated by the triangles, circles 
and squares respectively. 


t/T pa 1,2,3... (shaded regions in figure 3(b)), the capsule-wall distance h cw = h ini + /q at is 
at a local minimum and the background shear 7 reaches its maximum value. At this time, 
the background velocity 7 h cw , the product of the two, decreases because the phase difference 
is around 7 r. This phase difference does not exist in steady shear when there are no delays 
and the local background flow depends only on the capsule-wall distance h cw . 

The time-averaged value of the migration velocity U\ a t/a'y mgbX and of the slip velocity 
[4iip/a 7 max of the three capsules are indicated in figure |4] by horizontal lines. For all cases, 
the mean migration velocity is positive: the capsule undergoes a net migration away from 
the wall. However, the mean slip velocity is almost zero, namely, the capsule has no net 
motion in the streamwise direction. The capsule with the intermediate capillary number 
Ca = 0.15 has a higher migration velocity than its floppy (Ca = 0.0375) and stiff (Ca = 0.6) 
anmterparts. To further examine this dependence, the migration velocity Ci a t/a 7 m a x is 
therefore displayed versus the capillary number Ca in figure [5] for three different initial 
offsets. The mean velocity varies non-monotonically with Ca, the optimal capillary number 
Ca opt being around 0.1. This non-monotonic dependence of the migration velocity is in 
contrast to what is observed in steady shear, where the lift velocity increases monotonically 
with Ca^^. 


In steady flows, the particle lateral migration velocity is_a linear 


stress differences N± — N 2 as shown among others in Ref. 


60 and 


37 


unction of the normal 
. It is thus natural to 


investigate whether this relation holds in an oscillating shear as well. The normal stresses 
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Ni and N 2 are defined as, 


Ni = 


a. 


yy 


( 5 ) 


No = a 


yy 


<J, 


where cr denotes the extra stress due to the presence of the capsule in the volume V : 

on = -y J s Pi x jdS. 


( 6 ) 


We consider the nondimensional stress difference (Ah — N 2 ) //i0 v 7 max , with (f) v = 1/V the 
volume fraction; in the dilute limit, the normal stress difference AN = (N\ — N 2 ) //x0 v 7max 
represents the capsule-induced stress normalized by the viscous stress. 

The temporal evolution of the migration velocity Ci a t/a 7 max and the normal stress differ¬ 
ence AN are depicted in the top panel of figure [ 6 ] for capsules with Ca = 0.075, Ca = 0.12 
and Ca = 0.6. After the transient period, both quantities vary periodically and reach their 
local maxima with small phase delays, which indicates a considerable correlation between 
these two quantities. The period is around Tj 2, indicating that C/iat and AN change twice 
as faster as the oscillating shear. 

We next examine the local maxima of U\. At /a A f m ^ (marked by circles in the figure) and 
of AN (indicated by squares), always neglecting the peaks in the transient period. For the 
stiff capsule Ca = 0.075 (cf. figure 6 (a)), we identify five instants when f/i at /a 7 max and AN 


both attain local maxima. At times t \, £ 3 and £ 5 , the peaks of the two quantities have a 
similar magnitude with respect to their own scales, while at the time t 2 and £ 4 , their relative 
magnitude is significantly different. In other words, the ratio between the magnitude of the 
two local maxima varies with time, something which becomes more evident for the floppy 
capsule Ca = 0.6. In contrast, the two quantities U\ at /aj msix and AN have roughly the 
same magnitude every time they reach a local maxima for the capsule with Ca = 0 . 12 , 


the one displaying the largest lateral migration velocity (see figure 6 (b)): the ratio of the 
peak values remains almost constant in time. This implies that the correlation between the 
two quantities is stronger when Ca « 0.1 and decreases for larger and smaller values of the 
capillary number (Ca = 0.075 and 0.6 in the figure). This high correlation suggests that the 
normal stress difference contributes to a more pronounced lateral migration. 

To gain further insight, we show the time-averaged migration velocity Uia.t/aj max together 
with the time-averaged normal stress difference AN versus the capillary number Ca in 
figure [D the data reported pertain to three different values of the initial distance from the 
wall. Interestingly, these two quantities follow closely the same trend: a sharp increase with 
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FIG. 6. (Color online) Top panel: lateral migration velocity, Ci a t/a 7 max , (solid lines) and normal 
stress difference AN (dashed lines) versus the nondimensional time t/T, for three capillary numbers 


Ca = 0.075 (a) 0.12 (b) and 0.6 (c) Bottom panel: deformation index D versus t/T for the same 
capsule. The initial offset is h- in \/a = 2 and the frequency of the imposed shear ujj 7 max = 5/3. 
The local maxima of the velocity, stress and deformation are marked by the circles, squares and 
hexagons respectively; the shaded region indicates the period when they co-appear. 


Ca for stiff capsules, a maximum at Ca « 0.1 and then slower decrease for softer and softer 
membranes. 

The parameter study is continued by examining the dependence of the lateral migration 
velocity on the frequency of oscillation of the imposed shear, cu/ 7 max . Figure E] displays the 
time-averaged migration velocity, Ui a .t/a , j max , versus u>/ 7 max for the capsule with Ca = 0.3 
and the three offsets h ini /a = 1.5, 2 and 3. As cu/ 7 max < 2 , U\ a ,t/a'/ maiX decreases sharply 
with cn/ 7 max , almost linearly; for u>/ , y max > 2 , the migration velocity decreases more slowly, 
approaching asymptotically zero for large u>/ 7 max . The cases with cu/ 7 max = 0 correspond 
to a capsule in steady shear flow. Conversely, as co/'y mSLX oo, the flow oscillations are 
much faster than the relaxation time of the capsule and the capsule therefore does not have 
sufficient time to adjust to the flow before it changes direction: the deformation and the 
consequent lateral migration are hence negligible. 

The variations of the time-averaged migration velocity Ui a t/a A f max with the oscillation 
frequency are displayed in figure [9] versus the capillary number Ca for capsules with an 
initial offset h- m i/a = 2. The non-monotonic variation of U]- dt Ja A / max with Ca holds for all 
cases, but the capillary number Ca opt of maximum lateral velocity changes with the frequency 
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(a) (b) (c) 


FIG. 7. (Color online) Time-averaged value of the lateral migration velocity U\ at /a'y max and of the 
normal stress differences AN, versus the capillary number Ca. t7i a t/ a 7max and AN are indicated 
by the solid (resp. dashed) curve with circles (resp. diamonds), measured by the left (resp. right) 


scale. The frequency of the shear is tu/ 7 max = 5/3 and the results are shown for (a) h- m i/a = 1.5, 
|(b)| him/a = 2 and (c) h\ n \/a = 3. 



FIG. 8. (Color online) Time-averaged migration velocity E7i a t/ a 7maxj as a function of the frequency 
of the oscillating shear cu/ 7 max - The capillary number Ca of the capsule is 0.3 and three initial 
positions of the capsules are chosen, h- m \/a = 3 denoted by squares, h-^/a = 2 by circles and 
him/a = 1-5 by triangles. 

oj/7 ma x• The inset of the same figure shows that Ca opt varies almost linearly with the inverse 
of the shear frequency TWx/^- In fact, as the migration velocity, normal stress difference 
and capsule deformation vary at a frequency of roughly 2 cu/ 7 max (see figure. EJ), 2 cu/ 7 max 
can be considered as the effective frequency of the flow oscillation determining the capsule 
behavior. Accordingly, the ratio of the relaxation time of the capsule over the effective time 
scale of the flow is 2Cao;/7 max . 
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The inset indeed shows that 2 Ca opt a;/ 7 max ~ 0.5. In another words, lateral migration 
becomes strongest as the effective time scale of the flow is of the order of the relaxation time 
of the elastic membrane. 


It is worth noting that this scaling has similarities with that observed in many other 
situations where the flow dynamics is significantly influenced by the nonlinearity solely 
arising from the elasticity of the fluid or structure. A rotating helical slender body, a typical 
model micro-swimmer propelling like a cork-screw, attains the most efficient propulsion 
when the relaxation time of the polymeric fluid it is immersed in is of the order of the 
typical flow tim o 6 h 62 . Similarly, an elastic filamentr^^ or flapper— actuated in a viscous 
fluid for propulsion or pumping reaches maximum efficiency as the so-called sperm number—, 
the ratio of the time scale of the clastic structure over that of the flow, is of order one. We 
further note that the relation between the optimal capillary number and the frequency 
of the applied shear could be potentially used to aid the design of flow-assisted devices 
to sort deformable cells^I - — . Specific cells may be extracted from a dilute suspension in an 
oscillating Couette device if the frequency of the applied shear is tuned to match the capillary 
number of the targeted cells; these cells would in fact migrate with highest velocity. The 
frequency of the applied shear would have to depend on the deformability of the targeted 
cells, therefore the same device can be used to sort cells with different deformabilities by 
varying the frequency of operation. 


Furthermore, the relation 2 Ca opt c<;/ 7 max ~ 0.5 can also be implied theoretically in the 
limit of high-frequency oscillation, i.e. u;/ 7 max 1 , when the capsule undergoes small de¬ 
formation. We can assume the capsule reaches the peak deformation and migration velocity 
at the same capillary number Ca opt . In an unbounded oscillating shear, the leading order 
of peak deformation is^ D vea ,k = ^~r — and that of the equilibrium value D in the steady 
case is^£ D = ||Ca. Hence the capsule reaches the peak deformation and migration as 
2Ca op tCu/7max = 25 ~ 0.5. 


Finally, we report in figure flOl the time-averaged migration velocity Hi a t/a 7 max and the 
normal stress difference versus the capillary number Ca for two different oscillation frequen¬ 
cies, cu/d max = 0.5 and cu/ 7 max = 3. As shown in figured we can clearly identify a positive 
correlation between these two quantities. 
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FIG. 9. (Color online) The time-averaged lateral migration velocity Ui a .t/a A f max versus the capillary 
number Ca, for three shear frequencies, w/7 max = 0.5 (triangles), u/7 max = 5/3 (circles) and 
w/7max = 3 (squares), the initial offset h[ n [/a = 2. Ca op t indicates the capillary number of the 
capsule with the maximum migration velocity and the inset shows its dependence on 7 max /cu. 



FIG. 10. (Color online) The time-averaged value of the migration velocity Z7iat/°7 ma x and that of 
the normal stress differences AN versus the capillary number Ca. The initial offset is h- m \/a = 2: 
(a)j w/7 max = 0.5 and El W / 7max — 3. 


C. Deformation of the capsule 

A capsule typically evolves into a prolate when subject to an unbounded shear flow. 
Hence, its deformation is usually quantified by the so-called Taylor parameter D = (T max — 
-Zi min )/(Tmax + Tmin)^, where L max and L min are the length of the major and minor axis of the 
elliptical profile in the shear plane. In our case, the profile of the capsule on the shear plane 
is not a ellipse, since the symmetry is broken by the presence of the wall as also observed by 
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FIG. 11. (Color online) Temporal evolution of the deformation index D for capsules with initial 
offset hi n i/a = 2 and frequency of the shear u;/7 max = 5/3. Results for Ca = 0.03, Ca = 0.15 and 
Ca = 0.6 are indicated by the solid, dot-dashed and dashed lines respectively. The dotted curve 
denotes the oscillating background shear 7 max cos (ojt), arbitrarily scaled for a better visualization. 
The profiles of the capsules reaching the peak deformation are displayed on the shear plane. 

Nix et al.—. We thus introduce the deformation index, D = — L' min )/(L' max + L' min ), 

where L' mSLX is the maximum distance measured from the capsule surface to the its center of 
mass and L' min is the minimum distance. The deformation index D quantifies how far away 
the capsule is from its stress-free shape, a sphere in the current case. 

The time evolution of the deformation index D is illustrated in figure Qj] for three values 
of the capillary number, Ca = 0.03, 0.15 and 0.6, same initial offset, h ini /a = 2 , and shear 
frequency u: / 7 max = 5/3. The index D displays wake-like variations in time. An initial 
overshoot is observed for the soft capsules, Ca = 0.3 and 0.6, similar to that observed 
in the wall-bounded steady shear flow^-. After one or two periods, the amplitude of the 
oscillations becomes constant. The initial transient is more evident for the relatively soft 
capsules, Ca = 0.15 and 0.6; indeed, a higher capillary number implies a longer relaxation 
time, hence the capsule needs more time to adapt to the unsteady flow. We also note in the 
figure that the frequency of the periodically varying deformation index D is twice that of 
the background shear. The stiffest capsule, Ca = 0.03, displays the maximum deformation 
roughly after the local extrema of the shear and the minimum deformation as the flow 
reverses its direction. The stiff capsule therefore feels and responds to the variations of the 
flow fast enough to change its shape accordingly; following so closely the flow, its maximum 
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FIG. 12. (Color online) The temporal peak of deformation -D pea k versus the capillary number Ca, 
for an initial offset h in i/a = 1.5, 2 and 3. The shear frequency is u/ 7 max = 5/3. 

deformation appears very shortly after the instant of maximum shear. The opposite applies 
to softer capsules. 

The temporal evolution of the capsule deformation was also reported in the bottom panel 
of figure |6l The shaded regions clearly indicate that the local maxima of the migration ve¬ 
locity, normal stress difference and deformation index appear roughly at a same moment; 
moreover, the three quantities vary with twice the frequency of the background shear. Such 
a coincidence is not surprising; in fact, for a droplet and a capsule in steady shear, theoret¬ 
ical and numerical works have identified the relation between its migration velocity and 
deformation as U^t ~ D (a/h-mi) 2 - Our analysis further confirmed that even in the presence 
of flow oscillation, the deformation and migration velocity, which are time-dependent in this 
case, are correlated. 

To better understand the capsule deformation, we examine the peak deformation D pea k, 
the maximum value of D after the initial transient overshoot. Figure [T2] depicts D pea k as 
a function of the capillary number Ca for the same three initial heights considered above. 
-Dpeak increases almost linearly with the capillary number Ca for Ca < 0.1. This resembles 
the linear relation between the Taylor parameter D and the capsule capillary number in an 
unbounded steady shear flow. Further increasing the capillary number, -D pea k increases more 
slowly and reaches an asymptotic value close to 0.3 as Ca > 0.6. This is in contrast to what 
observed in an unbounded steady shear where the deformation D always increases with Ca 
although slowly at large Ca. 

In oscillating flows, the deformation not only depends on Ca but also on the frequency 
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(Color online) (a) the peak deformation D pea k as a function of the shear frequency 
for /ij n i/a = 1.5, 2 and 3 and capillary number Ca = 0.3. |(b)[ D pea k versus ( a/h- m i ) 3 for 
= 1, 5/3, 2 and 10/3. 


of the oscillations, cc/ 7 max . The time needed to reach the maximum deformation increases 
with Ca; for large values of Ca, the capsule fails to reach the maximum possible deformation 
that would occur in a steady shear before the flow changes direction if the oscillations occur 
fast enough. Considering T/2, the time during which the shear has constant direction, the 
deformation of the capsules is hence limited by T instead of by its deformability for large 
values of the shear oscillation frequency. The relation between the deformation and cu/ 7 max 
can also explain the non-monotonic dependence of the migration velocity with Ca. The 
lateral migration velocity, C/i at , of a capsule with A = 1 in wall-bounded shear is related to 
Smem, the VV component of the stresslet S mem induced by the elastic force on the membrane 
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certain value. In another words, the lateral migration of high-Ca capsules is hampered by 
the limiting deformation they attain. 

We further examine the dependence of D peak on the frequency of the shear cj/ 7 max and 
on the initial height h m i/a in figure fl3l where the membrane deformability is kept at Ca = 
0.3. The maximum deformation D peSL k decreases monotonically with uj/ 7 max as shown in 


figure 13(a) -D pea k increases as h ini /a decreases, i.e. larger confinement effects, and this is 
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FIG. 14. (Color online) Phase delay 0 of the capsule deformation with respect to the applied shear 
versus the capillary number Ca; the frequency of shear is uj/ 7 max = 5/3. Results are shown for an 
initial offset /i; n i/a = 1.5 (triangles), 2 (circles) and 3 (squares). 


more evident for smaller values of cu/ 7 max ; for large u>/7 max , the deformation barely varies 


with the initial height for the same mechanism explained above. We show in figure 13(b) that 
the maximum deformation D pea k is proportional to (a/h ini ) 3 . Note that in a wall-bounded 
steady shear flow, the same power law has been identified for the deformation of a droplet 
by Shapira and Haber— theoretically and for a capsule by Singh et a/.- 41 using numerical 
simulations. 


To quantify the delay between the capsule deformation and the shear oscillations, we 
define the time difference (t max D — t max s) scaled by the characteristic flow time as the phase 
lag 0 = (t m axD — tmaxs) /T. The phase delay 0 is negligible for the stiff capsule,whereas it is 
0 « 1/4 for the floppy capsule with Ca = 0.6; the latter deforms most when the flow changes 
direction (see the time histories in figure fTTT) . Note again that, regardless of the difference in 
the phase delay between capsules of different capillary numbers, the deformation oscillates 
at a frequency around 2 w/ 7 max . 

Finally, we examine in figure [14] the dependence of the phase delay 0 on the capillary 
number Ca. The phase delay 0 does not vary significantly with the initial offset h- m \/a of 
the capsules. 0 increases monotonically with Ca, an effect more pronounced when Ca is 
below 0.2. The delay 0 reaches an asymptotic value of approximately 0.25 as Ca continues 
to increase. Note that the phase delay 0 assumes a value between 0 (for purely viscous fluid) 
and 0.5 (for a linearly elastic solid). The asymptotic value of cs 0.25 clearly indicates that 
viscoelastic effect is an important feature in the flow of deformable cells. 


20 




IV. CONCLUSIONS 


Using an accelerated boundary integral method for the flow coupled with a global spectral 
method for the membrane, we study the near-wall dynamics of an initially spherical capsule 
with neo-Hookean membrane in an oscillating shear flow. We focus on the lateral migration 
and deformation of the capsule, and their dependence on the capillary number, the frequency 
of the background shear and the initial wall-capsule distance. 

The shape of capsule is asymmetric due to the presence of one wall and hence its center of 
mass and surface centroid are not equivalent, as discussed in Ref. [42. The time history of the 
lateral migration of the two centroids are clearly different, their time-averaged migrations 
being however similar. It is thus important to specify the definition of the centroids of 
asymmetric deformable particles when quantifying their migrations. 

The capsule reaches a quasi-periodic steady state after an initial transient period. During 
this state, the capsule follows a wiggling trajectory, moving away from and towards the 
wall periodically, with a net lateral migration away from the wall. For the capsule-wall 
distances and shear frequencies investigated here, the mean migration velocity varies non- 
monotonically with the capillary number Ca, reaching the maximum migration at an optimal 
capillary number denoted Ca opt . 

The optimal capillary number Ca opt is sensitive to the frequency of the shear, and is shown 
to scale linearly with the inverse of u;/ 7 max . Interestingly, it does not vary significantly with 
the capsule-wall distance. The maximum migration is observed when the effective oscillation 
period is of the order of the relaxation time of the elastic membrane. 

The capsule lateral migration velocity decreases monotonically with the frequency of the 
imposed shear; at relatively high oscillation rates, the capsule fails to adapt to the flow 
before it changes direction so that its lateral migration and surface deformation become 
negligible. 

The relation between the lateral migration velocity and the normal stress difference in¬ 
duced in the flow is explored. The maxima of these two quantities are closely correlated, 
with a weak time delay function of the capillary number. The relative magnitude of the cor¬ 
relations can vary significantly. As Ca ~ Ca op t, when the capsule has the highest migration 
velocity, the correlation is highest. The dependency of these two quantities on the capillary 
number Ca is investigated for different capsule-wall distances and oscillation frequencies of 
the imposed shear; in all cases, the two curves agree well with each and the peak values occur 
as Ca ~ Ca op t- This confirms the correlation previously observed in steady wall-bounded 
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shear flows 37 . It is worth pointing that such a relationship has also been discovered for 
a vesicle in two-dimensional unbounded Couette flow—; the vesicle migrates towards the 
center at a velocity linearly scaling with the normal stress difference. 

The deformation of the capsule exhibits a periodic variation approximately in phase with 
that of the migration velocity and of the normal stress difference, also at twice the frequency 
of the imposed shear. The maximum deformation increases linearly with Ca when Ca <0.1 
and more slowly as Ca >0.2. It reaches an asymptotic value when Ca is above a critical 
threshold Ca ~ 0.6; in this regime the deformation is limited by the time over which constant 
shear is applied and not by the membrane deformability. The peak deformation is found 
to scale with the capsule-wall distance (/i ini /a) 1//3 , as observed for capsules and droplets in 
near-wall steady shear—. We also discuss the phase delay between the capsule deformation 
and the background shear, another viscoelastic feature of the capsule dynamics. 

The non-monotonic dependence of the capsule migration velocity on its deformability 
could be potentially used to sort cells. Oscillating a suspension in a Couette device, cells 
with the capillary number exhibiting the largest lateral migration velocity may be extracted 
and isolated. By tuning the oscillation frequency, specific cells may therefore be targeted. 
The role of the oscillations on the rheological behavior of a suspension is a challenging 
extension of the present work. 
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